Origins of dihydrogen binding to metal-inserted porphyrins: 
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Q. ■ Abstract 

Using density functional theory calculations, we have investigated the interactions between hy- 
drogen molecules and metalloporphyrins. A metal atom such as Ca or Ti is introduced for in- 
^ '. 
"^ ' corporation in the central N4 cavity. Within local density approximation (generalized gradient 

c/3 '. 

D . approximation), we find that the average binding energy of H2 to the Ca atom is about 0.25 (0.1) 



C^ 



CJ 



eV/H2 up to four H2 molecules, whereas that to the Ti atom is about 0.6 (0.3) eV per H2 up to two 
H2. Our analysis of orbital hybridization between the inserted metal atom and molecular hydrogen 



'^ ■ shows that H2 binds weakly to Ca-porphyrin through a weak electric polarization in dihydrogen, 

O . but is strongly hybridized with Ti-porphyrin through the Kubas interaction. The presence of d 



orbitals in Ti may explain the difference in the interaction types. 
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I. INTRODUCTION 

Since hydrogen is a clean and recyclable energy carrier producing no pollution, there has 
been a great deal of effort to find good hydrogen storage materials. So far, much attention 
has been paid to low-dimensional nanostructures, since compact and light-element storage 
materials are required for mobile applications. In particular, carbon-based nanostructures 
such as carbon nanotubes, fullerenes and polymers were expected to be good candidates 
for hydrogen storage because of their large surface areas to adsorb hydrogen.-"- However, 
pure carbon-based structures have a drawback of insufficient binding energy^i^ near room 
temperature: their binding energies are lower than that for good hydrogen storage which 
should be 0.2 — 0.3 eV to operate at ambient conditions.-"— One the other hand, carbon- 
based nanostructures have been studied as promising catalysts for hydrogen release.—"— 
Berseth et alM showed that carbon nanostructures may be used as catalysts for hydrogena- 
tion/dehydrogenation of sodium aluminium hydride (sodium alanate; NaAlH2). 

As candidates having aforementioned optimal binding energies among nanostructures, 
metal-decorated materials such as metal-carbon based materials,-'^^"— metal-organic frame- 
works,—"— and metal-decorated polymer s^'22i'2£ were suggested to increase the potential for 
hydrogen storage. For cis- and trans-polyacetylene decorated with Ti atoms,— up to five 
and four hydrogen molecules binds to a Ti atom, respectively, and the optimal binding en- 
ergy of H2 molecules is around 0.2 — 0.3 eV. Because of their high cohesive energy (> 4 eV), 
however, transition metals (TMs) such as Ti and V have a strong tendency of clustering. 
Besides they can act as catalysts to destroy the carbon-based frameworks or to dissociate 
H2. Consequently, the hydrogen storage capability is significantly reduced. To prevent clus- 
tering of metal atoms, researchers have considered alkali-earth metals (AEMs) with lower 
cohesive energy (1 — 2 eV). Recently, Yoon et ai— showed that a coating of Ceo fullerenes 
with Ca could be a non-dissociative hydrogen strorage material, using ab initio calculations 
within generalized gradient approximation (GGA). On the other hand, it is known that alkai 
metal (AM) elements have relatively low cohesive energies (< 1.7 eV), which can also avoid 
clustering effectively.— Li-covered graphene were investigated to show that each adsorbed 
Li on graphene adsorbs up to four H2 molecules amounting to 12.8 wt%.— For alkali-metal- 
decorated organic molecules, Huang et al. found that Li-doped cases exhibit a higher storage 
capacity (> 10 wt%) than Na- and K-doped cases with the adsorption energy of ~0.1 — 



0.3 eV/H2.— Here, it is noted that AM elements are more reactive than AEM elements in 
ambient conditions. 

As a famous group in organic chemistry, porphyrin consists of four pyrrole-like subunits 
joined together by four methine (=CH-) groups. Four nitrogen atoms in a porphyrin molecule 
can be saturated by a single metallic ion instead of two H atoms. In nature, a heme contains 
a porphyrin coordinated to Fe, and a chlorophyll contains a porphyrin coordinated to Mg. 
Porphyrins are easy to synthesize and become important molecular building blocks that 
can be used to make several structures in materials chemistry.—"— Recently, many studies 
focus on the application of various properties of porphyrin.— i^ Previously, we also studied 
energetics of Ca-porphyrin where Ca is incorporated into the central N4 cavity of porphyrin.— 

In this paper, we present the different origins of dihydrogen binding to AEM- and TM- 
inserted porphyrin molecules on the basis of first-principles electronic structure calcula- 
tions. Ca and Ti were chosen as metal atoms inserted in the central N4 cavity of porphyrin 
(C20H12N4). First, we compared energetics of both Ca- and Ti-porphyrins. Then we inves- 
tigated interaction between dihydrogen and an AEM (Ca) or a TM (Ti) atom incorporated 
to porphyrin. As a result, we found that dihydrogen binds weakly to Ca atom through 
a weak polarization of dihydrogen, but interacts strongly with Ti atom with the so-called 
Kubas interaction. Finally, we discussed a possibility of porphyrin-based superstructures as 
hydrogen storage materials. 

II. CALCULATION METHODS 

We carried out first-principles calculations within local denstiy approximation (LDA) and 
generalized gradient approximation (GGA), using ultrasoft pseudopotentials implemented 
in the Vienna ab initio Simulation Package (VASP).— i^i The atomic positions were relaxed 
with residual forces smaller than 0.01 eV/A and the energy cutoff for a plane wave basis 
set was 400 eV. The vacuum region was about 10 A between porphyrin molecules, which 
was enough to avoid the image-image interaction. First, we optimzed the structure of 
metalloporphyrin such as Ca- and Ti-porphyrin molecules. The inserted metal (Ca or Ti) 
atom is located at the center of the cavity surounded by four N atoms. The Ca atom is 
placed above the porphyrin plane while Ti atom is located at the same plane to the z-axis 
in the metalloporphyrin. For the Lowdin charge analyses and the wave function plots, we 



employed the Quantum Espresso package. 

III. RESULTS AND DISCUSSION 

In a previous paper,— we briefly studied H2 adsorption on Ca-porphyrin to find that 
binding energy of H2 is about 0.25 eV per H2, using LDA up to four adsorbed hydrogen 
molecules. Here, we present the detailed results of binding energies and configurations for 
nB.2 (n = 1, 2, • ■ • , 6) on Ca-porphyrin along with the case of Ti-porphyrin. The binding 
energy Eh per H2 for n hydrogen molecules adsorbed on metalloporphyrin (M-porphyrin) is 
defined by 

Eb = (l/n)[£'(M-porphyrin) + nE(R2) - -E(nH2/M-porphyrin)], 

where £^ (M-porphyrin), E{B.2), and £'(nH2 /M-porphyrin) are the total energies of M- 
porphyrin, H2, and nH2 /M-porphyrin, respectively. The binding energies for Ca-porphyrin 
as well as Ti-porphyrin using both LDA and GGA are listed in Table I, while the respective 
optimized geometries are shown in Figs. [T] and [21 

For Ca-porphyrin, tilted configurations of H2 molecules are the most stable up to 4 
adsorbed hydrogen molecules. The GGA binding energies are quite small (about 0.1 eV/H2) 
compared to the LDA values (about 0.25 eV/H2 up to four H2). Note that CaCeo attracts 
up to 5H2 molecules with a uniform binding energy of ~0.2 eV/H2 for GGA and ~0.4 eV/H2 
for LDA.— Since the binding energies for the fifth and sixth molecules are relatively small, 
the H2 molecules are considered to be physisorbed on porphyrin. For our LDA results of 
Ca-porphyrin, the bond lengths of H2 molecules were slightly elongated to ~0.78 A from 
0.765 A (by ~2%), and the distances between H and Ca atoms were 2.35 — 2.52 A, which 
are very similar to the sum of atomic radii of H (0.53 A) and Ca (1.94 A). It means that 
there is no practical covalent bonding between H and Ca. 

Compared with results of Ca-porphyrin, the energetics of H2 adsorption on Ti-porphyrin 
was also studied. The LDA and GGA results for the H2 binding energy on Ti-porphyrin are 
summarized in Table I. Although the number of H2 adsorbed on porphyrin are only two, the 
binding energy is even larger than that for Ca-porphyrin. The LDA (GGA) binding energies 
per H2 are 0.63 eV (0.32 eV) up to two H2. 

Next, we considered a binding mechanism between adsorbed H2 molecules and the Ca 



atom. Though there are several energy levels aligned at the same positions in the PDOS 
of H2/Ca-porphyrin (Fig. |3]in Ref. l35l ) as well as 4H2/Ca-porphyrin, no significant overlap 
(i.e., orbital hybridization) was found between d orbitals of Ca and a orbitals of H2 in the 
wave functions. We calculated the charge difference, Ap, between the configrations before 
and after adsorption of a hydrogen molecule; 



^P = P„H 



2/M-pOrphyr: 



in ' 11H2 » M-porphyrin ' 
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From the Lowdin charge difference between total charge density of 4H2/Ca-porphyrin and 
the sum of those of Ca-porphyrin and 4H2, as shown in Fig. 3, we found that there is small 
electric polarization in each hydrogen molecule, which is associated with a weak interaction 
between H2 molecules and the Ca ion. According to our charge analysis, Ca donates about 
0.8 e to porphyrin, and becomes a positive ion (Ca"*"). Consequently, the tilting shape gives 
rise to the weak electric polarization in a H2 molecule under the electric field by the Ca ion. 
The average difference in charge in the two H atom in H2 is 0.015 e. It is hard to assign the 
number of electrons clearly to a specific atom using the Lowdin charge analysis as well as 
other charge analysis methods. Thus, we have to accept the value of charge with caution. 
This phenomenon was also reported in H2 binding to a Ca-hydroxyl group complex.— 

Then, we examined interaction between H2 and Ti-porphyrin. Similar to Ca-porphyrin, 
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FIG. 1: Optimized configurations of H2 molecules adsorbed on Ca-porphyrin. The upper and lower 
panels show top and side views of optimized structures, respectively. Yellow, blue, green, and white 
balls represent carbon, nitrogen, calcium, and hydrogen atoms, respectively. 
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FIG. 2: Optimized configurations of H2 molecules adsorbed on Ti-porphyrin. The upper and lower 
panels show top and side views of optimized structures, respectively. Yellow, blue, red-violet, and 
white balls represent carbon, nitrogen, titanium, and hydrogen atoms, respectively. 

we considered a charge difference between the total electron density of 2H2 /Ti-porphyrin 
and the sum of those of two isolated H2 and Ti-porphyrin. As shown in the charge difference 
plot in Fig. m an electric displacement around both Ti and two H2 occurs along the vertical 
line to the porphyrin. Such relatively large and complicated electric displacement is related 
to the Kubas interaction, which is roughly defined by a sort of interaction between d orbitals 
of metal atom and a (or a*) orbitals of H2.— In addition, it shows that charge redistribution 
takes place in the Ti orbitals (especially for d^^) because of adsorption of H2. For the case 
of Ti-porphyrin using LDA, the bond lengths of H2 were considerably elongated to 0.796 — 



TABLE I: The binding 
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0.810 A from 0.765 A (by ~ 4 — 6%), and the distances between H and Ti atoms were 1.9 
— 2.0 A, which are shorter than the sum of atomic radii of H (0.53 A) and Ti (1.76 A). It 
imphes that H2-Ti binding is stronger than H2-Ca binding. 

To understand interaction between H2 and Ti atom in more detail, we considered the 
electronic properties of the 2H2/Ti-porphyrin system such as the projected density of states 
(PDOS) and wave function characters. The PDOS of Ti, adsorbed H2 molecules and sur- 
rounding N atoms is shown in Fig. [5t^a), where some (apparent) hybridization are seen 
between d orbitals of Ti and s orbital of H (or N). After detailed investigations of wave 
function characters for the coincident energy levels, we obtained a hybridization diagram in 
Fig. Et^b), where we found the bonding and anti-bonding states between Ti "id orbitals and 
(J (or 0"*) orbitals of H2. Each wave function character labeled ^a'-^d^ corresponds to each 
denoted electronic states having the same label in Fig. [5t^a).— 

Three hybridizations labeled ^a''-^d are clearly noticeable. The first one denoted by 'a' is 
a hybridized bonding state of the empty d-^2 orbital of Ti with filled a orbitals of H2, where 
the filled a orbital donates electrons to empty (i^2 orbital. The second one denoted by '6' 
just below the Fermi energy is another hybridized bonding state of the filled d orbitals {dzx, 
dzy) of Ti with empty a* orbitals of H2, where the filled d orbitals {dzx-, dzy) back-donate 
electrons to empty a* orbitals. This bonding state has a 'twin' anti-bonding state denoted 
by '6' which is located just above the Fermi level: Due to surrounding environments around 
Ti-H2, these two states {b and b') would be split from an original bonding state hybridized 




FIG. 3: Isosurface plot of charge difference between the total electron density of 4H2/Ca-porphyrin 
and the sum of those of isolated 4H2 and Ca-porphyrin. Charge accumulation and depletion are 
represented by red and blue, respectively. The isosurface levels are ±0.0025 electrons/bohr^. A 
slight electric displacement of hydrogen molecules is shown. 




FIG. 4: Isosurface plot of charge difference between the total electron density of 2H2/Ti-porphyrin 
and the sum of those of isolated 2H2 and Ti-porphyrin. Charge accumulation and depletion are 
represented by red and blue, respectively. The isosurface levels are ±0.0025 electrons/bohr^. Com- 
plicated and relatively large electric polarization is shown around Ti and H2 molecules. 

from the filled d orbitals {dzx, dzy) of Ti with empty a* orbitals of H2. On the other hand, we 
have a distinct anti-bonding state (denoted by 'c') obtained from the (1^2 — a hybridization. 
Another anti-bonding state between d orbitals [d^^^ dzy) and a* of H2 is not clearly shown. 
As a result, we conclude that H2 binds strongly to Ti atom through the Kubas interaction. It 
should be compared with a weak interaction between H2 and Ca through slight polarization 
of hydrogen molecules. A difference in the interaction types between Ti- and Ca-porphyrin 
molecules may be attributed to the presence of d orbital in Ti. 

Finally, we turned our interest to a possibility of metalloporphyrin as a hydrogen stor- 
age material by calculating the hydrogen storage capacity. When hydrogen molecules are 
directly adsorbed on the Ca atom, the storage capacity is estimated to about 2.3 wt% for 
4H2/Ca-porphyrin. On the other hand, the hydrogen storage capacity of Ti-porphyrin is 
about 1.1 wt% for 2H2 /Ti-porphyrin. Unfortunately, these values are much smaller than 
the requirement for practical applications of hydrogen storage. Whereas model systems 
of other research groups are carbon-based structures with metal adatom, our model has 
strong covalent bonding between the Ca (or Ti) atom and four nitrogen atoms in porphyrin. 
Thus, there are not so reactive electronic orbitals to interact with H2 in metalloporphyrin. 
Consequently, the storage capacity of metal-incorporated porphyrin is not as high as those 
of metal-decorated carbon nanostrcutures. If metal organic frameworks (MOFs) are con- 
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structed using metal-incoporated porphyrins,— the storage capacity could increase even in 
ambient conditions. Therefore, further investigations are needed for metal-incorporated 
porphyrin molecules to seek high-capacity H2 storage materials. At present, gravimetric 
storage capacities of MOFs are much smaller than 1 wt% at ambient pressures and room 



temperature." 



42^ 



IV. CONCLUSION 

In summary, we have investigated non-dissociative molecular hydrogen binding onto M- 
inserted (M = Ca, Ti) porphyrin, and have compared the differences in the origins of dihy- 
drogen binding to Ca- and Ti-porphyrin, using density functional theory calculations. By 
successively increasing the number of adsorbed hydrogen molecules on a Ca-porphyrin or 
Ti-porphyrin, we obtained the optimized structures and the H2 binding energies. According 
to our LDA (GGA) calculations, the binding energy of H2 to the Ca atom is about 0.25 




FIG. 5: (a) PDOS and (b) hybridization diagram for 2H2 /Ti-porphyrin. The dotted line in the 
right panel represents the Fermi level. 



(0.1) eV/H2 up to four H2 molecules while that to the Ti atom is about 0.6 (0.3) eV/H2 up 
to two H2. Our results demonstrate that a calcium atom incorporated into porphyrin with 
H2 molecules forms a dihydrogen complex mainly through a weak polarization of H2. In 
contrast, Ti-porphyrin generates a Kubas-type dihydrogen complex with H2. The presence 
of d orbital in the Ti element may explain the difference in the interaction types. 
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